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Abstract 

We outline a non-perturbative procedure for calculating the total photoionization cross-section 
of two-electron atomic systems. The procedure is based on the Floquet-Fourie representation of the 
solution of the time-dependent Schrodinger equation. The Floquet-Fourie ansatz produces a set of 
equations which is recast into a generalized eigenvalue problem by means of the complex rotation 
procedure. With the use of the Hylleraas-type basis functions, the total photoionization cross- 
sections are obtained within the accuracy of a fraction of a percent. The total photoionization cross- 
sections for neutral helium are in good agreement with the convergent close-coupling calculations 
of Kheifets and Bray [Phys. Rev. A 58, 4501 (1999)] but deviate notably from the experimental 
data of Samson et al. [J. Phys. B 27 887 (1994)]. 
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I. INTRODUCTION. 



Photoionization of two-electron atoms has been studied theoretically by different authors 
starting from the pioneering paper by Wheeler Review of early literature on this subject 
can be found in Ref. [^J. Subsequently, a large number of computations of helium photoion- 
ization cross-sections was reported 3, J 0| These calculations produced a collection 
of results varying typically by 5% from each other. On the experimental side, the bench- 
mark set of data was reported by Samson et al. 3] who measured the total photoionization 
cross-section of He in the photon energy range from the threshold to 120 eV. Agreement 
between the theoretical and experimental data was within the same margin of 5%. In the 
following years, the theoretical interest shifted towards calculation of differential character- 
istics of the photoionization process and to studies of double photoionization. Hej^e^several 

con- 



approaches have been advocated including the many-body perturbation theory 
vergent close-coupling, method P, [llj], time-dependent close-coupling method [hi 

HQ 



3Q, 



i?-matrix approach |U,ll£|], and methods based on the computation of the dipole response 
function |16| or 5-spline implementations of the exterior complex scaling jl7|. 

Due to this shift of focus, there have been no further attempt to produce a consistent 
set of photoionization cross-sections of He with an accuracy better than several percent. 
Yan et al. ^| combined measurements of Samson et al. Q] at low energies and theoretical 
calculations at high energies to construct a set of photoionization cross sections of He that 
should be reliable at all energies. However, there was no consistency check applied to the 
experimental data. In the meantime, accurate helium photoionization cross-sections would 
be highly desirable due to importance of He in astrophysics and its use as a standard gas in 
determination of the photoionizaiton cross-sections of other atomic and molecular species. 

In the present paper, we develop the complex rotation method (CRM) for highly accurate 
calculations of the total photoionization cross-section of two-electron atomic targets. One 
way of calculating the photoionization cross-section is to combine the CRM technique with 
the perturbation theory with respect to interaction of the atom with the electromagnetic 
field. In such a perturbation theory, the CRM provides the basis of the field-free atomic 
states. It was demonstrated by Johnson and Reinhardt Q] that relying on the spectrum of 



the CRM eigenvalues, one can construct a representation of the complete Green's function of 
the atom. This, in turn, allows to write down a convenient representation for the projection 
operator corresponding to the continuous spectrum of the atom [20]. Using this projection 
operator, one can compute probabilities of transitions into continuum under the action 
of some perturbation, in particular, the interaction of the atom with the electromagnetic 
field. Calculations of total photoionization cross-sections of the helium atom based on this 
technique have been reported in Refs. 2jJ, l22j. Similar ideas were also used to determine 
static and dynamic polarizabilities of helium |23l. \24 \ . 

In the present work, we use the CRM procedure in a somewhat different, non-perturbative 
way by applying it to the whole system the atom plus the electromagnetic field. Thus, we 
are capable of going beyond the perturbation theory and considering very strong fields. 
In this respect, the present technique has certain features in common with Refs. 25, 26], 
where an approach based on the configuration interaction procedure was developed for atoms 
with more than one electron. We cast our formalism using the language of square-integrable 
functions with a finite norm. This approach becomes feasible in the so-called Floquet-Fourie 
representation of the time-dependent Schrodinger equation. 

Another key ingredient of the present work is the Hylleraas basis functions which have 
long been used in various variational-type calculations. An excellent review of applications 
of the Hylleraas basis to calculations of energies of two-electron atoms is given by Drake 
j^ . A well-known trademark of the Hylleraas basis set is a very high accuracy of the 
atomic energies. In the present paper, we show that the same high accuracy which is 
achieved for field-free atomic states can also be attained when the atom is subjected to a 
monochromatic electromagnetic field. In particular, the total photoionization cross-sections 
can be calculated with an accuracy of the order of a fraction of a percent. 

Thus generated cross-sections were compared with the experimental results of Samson 
et al. 3]. We discovered a systematic deviation from the experiment, especially in the 
region close to double ionization threshold at the photon energies of ~80 eV. This deviation 
was confirmed by comparison with earlier results produced by the convergent close-coupling 
(CCC) method ^ . 
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The true potential of the present method is realized in the strong field regime where 
the perturbation theory fails. As demonstrated below, the Floquet-Fourie-Hylleraas ansatz 
produces very accurate results in this regime as well. 



II. THEORY 



A. General Theory. 



The non-relativistic Hamiltonian of the helium atom in the presence of the external 
monochromatic linearly-polarized electromagnetic field can be written as : 



where, T is a kinetic term: 



U potential energy term: 



H = f + U + V, (1) 

f = ^- + ^-, (2) 

2 2' v ' 

2 2 1 

U = ---- + -^-, (3) 
7*1 r 2 |t*i — r 2 \ 

and V describes interaction of atom and the field. In the length gauge (which will be used 
in the present paper), this operator assumes the form: 

V = F D cos cut, (4) 

with D = r\ + r 2 . Unless stated otherwise, the atomic units are used throughout the paper. 

We write the solution of the time-dependent Schrodinger equation (TDSE) using the 
Floquet-Fourie ansatz 

tf(t) = e - iEt Y,u n e- inwt ■ (5) 

n 

By substituting this expression into the TDSE and equating coefficients with e - lEt - imwt ; we 
obtain a chain of coupled equations for the Floquet-Fourie coefficients u n \ 

F D 

(E -T -U + nu)u n = — - — (u n ^i + u n+1 ) , n = 0, ±1 . . . , (6) 
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We solve this set of equations with the help of the complex rotation procedure 

Formally, the CMR can be described as a complex transformation of radial variables 



r,i — > r^e 10 , where 9 is the rotation angle, the sole parameter defining the transformation. 



Under this transformation, the chain of equations © is converted into 

(E - f e- 2ie - Ue- ie + nu)u n = —^ e ie {u n - X + u n+l ) , n = 0, ±1 . . . , (7) 



According to the general theory of CRM the set of equations (J2J) can be 

solved by means of variational techniques if the rotation angle 9 is properly chosen. 

We introduce a basis set of square integrable functions \n, k) where the index n refers 
to the number of the Floquet block and the index k denotes a particular L 2 function in 
the subspace of the n-th block so that u n = ^c„fc|n, k). With these notations, the set of 

k 

Eqs.flJJ) can be rewritten in a matrix form as: 

F- D 

J2(ni,h\E + riw-Te- 2ie -Ue- i0 \n,k)c nk = ]T (m, h\— — e ie |n2&2>Cn 2 fc, (8) 

k 7i2=n±l,k " 

Notations can be further simplified by introducing obvious shorthands: 

pie 

((E + nuj)R: k iki -T^ ki e- 2W -U^ kl e~ w )c nk = £ V^f-c^ (9) 

n2=n±l Z 

where it is understood that summation is carried over the repeated k-index. Here V = 
F D, and R, T and U stand for the overlap, kinetic energy and potential energy matrices, 
respectively. 

One could say here a few words about the choice of the basis allowed by the structure 
of the system Q. Suppose first, that in each of the subspaces corresponding to different 
Floquet blocks we chose some compete set of functions, such that for any u n in Eq.(0) we 
had: u n = c nm \n, m). Let the set of |n,m)'s be the same for all Floquet subspaces. 
Then, if we have retained iV Floquet blocks in the system (J7J) and keep M terms in the 
expansion for each u n in Eq. (f7jl we have altogether NM unknowns c nm in the system (J2J). 
To get a correctly posed eigenvalue problem, we should have the same number of equations. 
This number is provided by projecting each of the equations ((Zj) on one of the \n, m)'s 
with m = 1 . . . M. This way of reducing the set of equations (JTj) to the form of matrix 
eigenvalue problem is correct, but too general for our purposes. It can be seen, that one can 
considerably diminish the resulting dimension of the matrix eigenvalue problem by using 
certain symmetry properties of the system Eq.fJZJ). It is easy to see, that this system allows 
the following class of solutions: m„'s with even n are of even parity, while u n 's with odd 



n are of odd parity. Parity here is understood with respect to the spatial inversion. Of 
course, there is a class of solutions with the opposite property: w n 's with even n are of odd 
parity, while u n 's with odd n are of even parity. The solution we are looking for (which is 
to describe behavior of the even l S e state of helium) evidently belongs to the first class. We 
can therefore, choose the basis set as follows. 

Instead of choosing the same set |n, m) for each Floquet block, we choose two sets: a 
set |n even ,m), consisting of basis finctions of even parity, is used as a basis to represent 
UnS with even n's. Another set |n odd ,m), composed of odd parity functions is used as a 
basis to represent w n 's with odd n's. Suppose that in the expansions of w n 's with even n's 
we retain M even terms, and in the expansions of m„'s with odd n's - M odd terms. Let the 
number of Floquet blocks with even and odd n's be respectively N even and N odd . Than 
we have iV even M even + N odd M odd unknown coefficients c nm . We obtain the same number 
of equations by projecting equations (J7|) on |n even ,m),m = 1. . . M even for even n and on 
l^odd, m),m — 1 . . . M odd for odd n. Projection of equations with even n on the |n odd , m) 
and of equations with odd n on the |n even ,m) gives identically zero and does not add new 
equations. More details about the basis functions Ineven,"^) and |n odd ,m) is given below. 

According to the general theory of CRM, some of the energy values (generally complex) 
for which system (JHJ) has a solution are related to the position and width of the resonance 
state via E = E r — iT/2, where E r is position of the resonance and T its width. This leads 
one to solving a generalized eigenvalue problem. Effectiveness of finding eigenvalues of such 
a problem depends crucially on the choice of the basis used to represent the matrices in 
Eq.flHJ). 

B. Basis set. 

The basis set used in the present paper was constructed from the Hylleraas type functions: 
9n,,n 2 ,N{rx,r 2 ) = rT r? |n - r 2 \ N e- ar ^ br2 \h(l)l 2 (2)L), (10) 
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where a,b are some constants (to be specified below), n%,n2,N are integers and the angular 
part 

|/i(l)/ 2 (2)L) = J2 fi^W^W, (11) 

represents two angular momenta Z l5 Z 2 coupled to a state with a total angular momentum 
L. The basis functions ()10|) must also be properly symmetrized with respect to exchange of 
the electron coordinates. When choosing parameters in Eq. (|l(Jj). we followed the following 



empirical rules [2JJ,|36j. All the basis functions with the parameters satisfying: 

n x + n 2 + N < A max (12) 

were included in the calculation. The parameter A max determines the overall size of the 
basis. There is also a semiempirical rule for choosing angular momenta I1J2 in the Eq.(|l(Jj). 
Thus, for states of the natural parity Zi are best chosen so that h + h = L. Both these 
criteria help to avoid the numerical problems due to near-degeneracy of the basis set when 
its dimension becomes large. 



III. NUMERICAL RESULTS 
A. Field-free case 

In the present work, our main goal is to obtain accurate photoionization cross-sections 
from the ground state of neutral helium for not very large electromagnetic field intensities. 
Accordingly, our main interest is focused on the states of S and P symmetries. Therefore, 
our first goal is to choose such a basis that solution of the eigenvalue problem yields 
accurate energies for the ground 1 S and first excited 1 P° state of the helium atom in the 
absence of the field. 

This goal was achieved as follows. We chose parameters iV max = 18, a = b = 2 for the 
^-states and N max = 13, a,b = 1,2 for the P-states. The reason for enlarging the basis 
set for the excited P-states is that the electrons in such states are generally on different 
distances from the nucleus. This choice combined with restriction on angular momenta ()12|) 
resulted in Ns = 372 basis functions for the ^-states and Np = 660 basis functions for the 
P-states. 
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The next step was to solve the generalized eigenvalue problem for the field-free case. In 
Eq. Q we put F = 0, u = 0, and limited ourselves to the blocks with n — 0, n = ±1, the 
n = block being composed of the states of 1 S e symmetry, and n — ±1 blocks composed of 
the states of l P° symmetry. All the numerical results reported below were obtained using 
the quadruple precision arithmetics. 

We note, that in the presence of the weak electromagnetic field account of the blocks 
with n — ±1 corresponds to absorption and emission of one photon. We shall use this fact 
below to extract the photoionization cross-section from our calculation. For the moment, 
we are concerned with testing the accuracy of our basis. Diagonalization of the eigenvalue 
problem with F = 0, uj = in the basis described above produced the following results 
for the complex energies: E = —2.903724384 + i 1.3 x 10~ 8 (the ground state) and E = 
—2.123843094 + i 7.6 x 10 -9 (ls2p 1 P° state). A small imaginary part which, in the absence 
of the field, should be zero could be taken as an indication of an accuracy of our basis set. 
Either this criteria or a direct comparison with the well-known results of highly accurate 
calculations |27J shows that we have achieved an accuracy of the order of 10 -8 a.u. This 
accuracy, as will be demonstrated below, is sufficient to obtain the photoionization cross- 
sections with at least three significant figures. 



B. Total photoionization cross sections 

To calculate the total photoionization cross sections we adopted the following strategy. 
The eigenvalue problem was solved with the Floquet blocks n = 0, ±1 retained, the 
composition of each block was the same as described above for the field-free case. Diagonal- 
ization of the eigenvalue problem produced energy shift and total width for the ground 
state. By definition, the photoionization cross-section from this state is related to the total 
width T via 

a = lim 87caTu/F 2 , (13) 

where F is field strength, uj its frequency, a is the fine structure constant. We need therefore 
to extract from our calculation the coefficient with F 2 in the asymptotic law defining the 
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TABLE I: Results for the ground state eigenvalue of problem @ as functions of parameters N x 



in Eq.JHJ), u = 80 eV, F = 0.1 a.u. 



N s 

* max 


N p 

A ' max 


Total dimension of the 
eigenvalue problem @ 


Re#(a.u.) 


T (a.u.) 


17 


11 


1300 


-2.90307660 


0.000487738 


18 


12 


1692 


-2.90307661 


0.000487698 


19 


13 


2204 


-2.90307659 


0.000487689 



weak-field behavior of the width: 

T(F) = T F 2 + TiF 3 + ... (14) 

To implement this strategy, we need an extrapolation procedure since the calculation based 
on the system (JHJ) is performed for a non-zero field strength. Although finite, this field 
strength should not be too small to compute T with sufficient accuracy. 

The issue of accuracy can be addressed as usual in variational-type calculations, by merely 
increasing the basis size and verifying that the results do not change appreciably. Such a test 
was performed for a photon energy ui = 80 eV and a field strength F = 0.1 a.u. by varying 
the parameter N max in Eq.(|12jl for the S and P states. The diagonalization of the problem 
© was performed with the Floquet blocks n — 0, ±1 retained. All the remaining details of 
the basis (nonlinear parameters etc.) were the same as in the field-free case reported above. 
The calculation was performed for the value of the rotation angle 9 = 0.3. 

The test results are presented in Table HI One can observe that, just as in the field- 
free case, the accuracy is on the level of 10~ 8 a.u., which implies that T has at least four 
significant digits in this interval of field strengths. 

The issue of the stability of the results with respect to the number of the Floquet blocks 
included in diagonalization of (jUJ) is addressed in the next section where we consider effects of 
going beyond the first order perturbation theory. We shall say in advance that including the 
Floquet blocks with n = ±2 in diagonalization of (jHI) does not alter the numerical accuracy 
appreciably. 
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TABLE II: Extrapolation of the T's to the zero-field limit. 



u (eV) 


F = 0.07 a.u. 


F = 0.1 a.u. 


T/F 2 (a.u.) 

F = 0.13 a.u. 


F = (Extrapolation) 


40 


0.4208622 


0.4201601 


0.4192063 


0.4215215 


80 


0.0488002 


0.0487698 


0.0487239 


0.0488112 


85 


0.0392854 


0.0392618 


0.0392330 


0.0393202 


91 


0.0306858 


0.0306720 


0.0306524 


0.0306961 


95 


0.0262180 


0.0262082 


0.0261936 


0.0262224 


111 


0.0147116 


0.0147084 


0.0147033 


0.0147116 


205 


0.0013719 


0.0013726 


0.0013729 


0.0013687 



As to the extrapolation procedure needed to extract the coefficient T in Eq. (fT4"|) . we 
chose a scheme based on the three-point Lagrange formula. For each frequency reported 
below, we performed calculations for the field strengths F = 0.07, 0.1, 0.13a.u. We also used 
a mid size basis set with A^ ax = 18, A^ ax = 12, Floquet blocks with n = 0, ±1, all other 
details of the basis are the same as in the field-free case above. Results of this calculation 
and extrapolation are shown in Table ITT1 

Using an estimate for the remainder of the series (fl"4"j) . it is a simple matter to verify that 
for the field strengths considered the possible relative error introduced by the extrapolation 
of T/F 2 is of the order of 0.1%. Hence, at least three digits in our result for the extrapolated 
ratio T/F 2 and the cross-sections reported below must be reliable. This level of accuracy can 
easily be improved by merely going to extrapolation schemes of higher order and computing 
T for more field values. 

In Table InTl we present our results for the cross-sections based on formula (|13|) in which we 
fed the extrapolated ratios from the last column of Table ILT1 Along with our data, we present 
the benchmark experimental results of Samson et al. Qlas well as earlier theoretical results 
from Ref . 2s| . The experimental setup of Samson et al. p| was such that the measured cross- 
section was summed over all final states of the remaining ion including the doubly ionized 
states. It is exactly the cross-section that is calculated presently and therefore comparison 
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TABLE III: Comparison of the present results and other theoretical and experimental data for 
the total photoionization cross section (in Mb). 



u (eV) 


Present 


L 


V 


ccc 

A 


Average 


Experiment 
Samson et al. [7] 


40 


3.1822 


3.188 


3.178 


3.247 


3.2043 


3.16 


80 


0.7369 


0.7432 


0.7403 


0.7366 


0.7400 


0.693 


85 


0.6308 


0.6364 


0.6327 


0.6294 


0.6328 


0.595 


91 


0.5272 


0.5333 


0.5284 


0.5248 


0.5288 


0.502 


95 


0.4701 


0.4765 


0.4717 


0.4689 


0.4723 


0.450 


111 


0.3082 


0.3097 


0.3089 


0.3081 


0.3089 


0.300 


205 


0.0529 


0.0533 


0.0534 


0.0531 


0.0533 


0.0510 



between the theory and experiment should be straightforward. 

n 

Our theoretical results agree with the data of Samson et al. p| within the postulated 
experimental accuracy of few percents. The strongest deviation is for uo = 80 eV where the 
difference between the present result and the experimental value is 6%. This is deviation is 
clearly seen in the Figure where we plot the present Floquet calculation along with the CCC 
calculation in three gauges of the lectromagnetic interaction and the experiment. Agreement 
between the present calculation and that of the CCC is much better, difference of the results 
of two approaches not exceeding 1%. The accuracy of the CCC result is hard to estimate 
directly as this method relies on the numerical solution of a set of close-coupling equations. 
The only implicit indication is the difference between the cross-sections calculated in the 
three gauges of the electromagnetic interaction, the length (L), velocity (V) and acceleration 
(A). This difference is typically 1-2%. Thus, the deviation of the present calculation with 
the CCC is more likely to be the problem of the latter as the former is believed to be much 
more accurate. 

As a by-product of the calculation described above, we also obtained the shift of the 
ground state of helium due to the interaction of atom with the linearly polarized monochro- 
matic field (Table ITVjl. 
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Floquet • 




80 85 90 95 100 105 110 115 120 

Photon energy (eV) 

FIG. 1: The total photoionization cross-section (sum of the single a + and double a 2+ cross- 
sections) as a function of the photon energy. The present calculation for selected photon energies 
is denoted by dotes. The CCC calculation in the three gauges of the electromagnetic interaction 
(Length, Velocity and Acceleration ) is exhibited by different colors / line styles. The experiment 
of Samson et al. |7| is presented by a dotted line. 

C. Extended calculation 

We now turn to extended calculations with inclusion of a larger number of the Floquet 
blocks n = 0, ±1, ±2 in Eq. Q. The aim of these calculations is two-fold. First, we shall 
confirm the stated accuracy of the present weak field results which is not effected by the 
number of the Floquet blocks retained in the calculation. Second, we report some preliminary 
results concerning behavior of the widths parameter in stronger fields where inclusion of a 
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TABLE IV: Real part of the energy of the ground state of helium in the presence of the linearly 
polarized monochromatic field. 



to (eV) 


F = 0.07 a.u. 


ReE (a.u.) 
F = 0.1 a.u. 


F = 0.13 a.u. 


40 


-2.90281954 


-2.90187690 


-2.90060016 


80 


-2.90340686 


-2.90307658 


-2.90263014 


85 


-2.90344158 


-2.90314741 


-2.90274973 


91 


-2.90347684 


-2.90321933 


-2.90287117 


95 


-2.90349708 


-2.90326062 


-2.90294089 


111 


-2.90355845 


-2.90338582 


-2.90315235 


205 


-2.90367802 


-2.90362977 


-2.90356450 



larger number of the Floquet blocks becomes essential due to a non-perturbative nature of 
the processes involved. 

The basis for the extended calculations was constructed as follows. As we discussed above 
the basis subset spanning each Floquet block in the system Q can be chosen to consist of 
the functions of a given parity, two adjacent blocks having opposite parities. Thus, in the 
low-field calculations described above, the block n = was composed of even basis functions 
while two blocks with n = ±1 contained odd basis functions. Inclusion of the blocks with 
n = ±2 is, therefore, equivalent to adding more even basis functions. We did it in the 
following way. In addition to the 1 S' e states we previously had in the n = block, the states 
of the symmetries 1 D e and 1 P e were included in the calculation. Thus the blocks with n = 
and n = ±2 had the following composition: A" max = 18 for the 1 S' e -basis functions, N max = 8 
for the 1 P e and 1 D e -basis functions. As before, the blocks with n = ±1 were composed 
of basis functions of 1 P°-symmetry with N max = 13. Thus, the basis set is considerably 
enlarged comparing to the one used in the previous section. With this choice of parameters 
N max , the overall dimension of the eigenvalue problem (jUJ) was 2676. Results produced for 
the ground state of He by diagonising this eigenvalue problem are shown in Table 

Comparison of the results given in Table IVl supports the assertion we made in the previous 
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TABLE V: Results of the calculation with Floquet blocks n = 0, ±1, ±2 included in the system 





0J 


= 111 eV 




= 205 eV 


F (a.u.) 


ReE (a.u.) 


T/F 2 (a.u.) 


ReE (a.u.) 


T/F 2 (a.u.) 


0.10 


-2.90338569 


0.014714 


-2.90362976 


0.0013734 


0.13 


-2.90315198 


0.014715 


-2.90356447 


0.0013739 


0.20 


-2.90236955 


0.014715 


-2.90334589 


0.0013743 


0.50 


-2.89525524 


0.014706 


-2.90135842 


0.0013745 


1.0 


-2.86985102 


0.014665 


-2.89426254 


0.0013701 



section as to the accuracy of our results for the widths. As one can see, for the field strengths 
F 0.1 a.u., inclusion of the additional Floquet blocks and basis states of symmetries other 
than S and P produces relative variations in the widths on the order of 0.01 percent. 
This means that for such field values we are still within the domain of the validity of 
the perturbation expansion. For the frequencies presented in the Table the domain of the 
perturbation theory actually extends quite far in the region of large field strengths. As one 
can see from the Table (0), the ratio T/F 2 starts changing in a more or less appreciable 
manner only for field strengths as large as F ~ 1 a.u. 



IV. CONCLUSION 

We performed a calculation of the total photoionization cross-sections from the ground 
state of helium. We employed a theoretical procedure based on the Floquet-Fourie repre- 
sentation of the solutions of the TDSE describing the helium atom in the presence of the 
linearly polarized monochromatic electromagnetic field. The resulting set of Floquet equa- 
tions was cast into a generalized eigenvalue problem by the complex rotation method. Our 
approach is essentially non-perturbative. This is in contrast with other works (20, 22 1 where 
the CRM was used to produce an accurate description of the field-free helium atom, thus 
giving the ground for application of the perturbation theory. In our approach, we do not rely 
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on any perturbation expansion to describe interaction of the atom and the electromagnetic 
field. This interaction is included into the theory from the beginning. We would like to 
emphasize the accuracy of the present results for the photoionization cross-sections which, 
we believe, is on the level of a fraction of a percent. Although only few selected photon 
energies were reported in the paper, far wider and denser energy grid was covered by the 
present calculation. These results might serve as an accurate database and find their use in 
various astrophysics and atomic physics applications. The authors shall gladly communicate 
these data on request. 
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